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Abstract. Simulations based on experimental data obtained from multifragmenting quasi- 
fused nuclei produced in central 129 Xe + nat Sn collisions have been used to deduce event by event 
freeze-out properties on the thermal excitation energy range 4-12 AMeV. From these properties 
and temperatures deduced from proton transverse momentum fluctuations constrained caloric 
curves have been built. At constant average volumes caloric curves exhibit a monotonous 
behaviour whereas for constrained pressures a backbending is observed. Such results support 
the existence of a first order phase transition for hot nuclei. 



1. Introduction 

One of the most important challenges of heavy-ion collisions at intermediate energies is the 
identification and characterization of the nuclear liquid-gas phase transition for hot nuclei, 
which was earlier theoretically predicted for nuclear matter pQ. In the last fifteen years a big 
effort to accumulate experimental indications of the phase transition has been made. Statistical 
mechanics for finite systems appeared as a key issue to progress by proposing new first-order 
phase transition signatures related to thermodynamic anomalies like negative microcanonical 
heat capacity and bimodality of an order parameter [2 [31 Hj. Correlated temperature and 
excitation energy measurements, commonly termed caloric curves, were the first studied possible 
signatures of phase transition. However in spite of the observation of a plateau in some caloric 
curves, no decisive conclusion could be extracted [H [5j El [7] . The reason is that experimentally 
it is not possible to explore the caloric curves at constant pressure or constant average volume, 
which is required for an unambiguous phase transition signature. Indeed, theoretical studies 
show that if many different caloric curves can be generated depending on the path followed 
in the thermodynamical landscape, constrained caloric curves exhibit different behaviours in 



presence of a first order phase transition: a monotonous evolution at constant average volume 
and a back bending of curves at constrained pressures [9]. With the help of a simulation able 
to correctly reproduce most of the experimental observables measured for hot nuclei formed 
in central collisions (quasi-fused systems, QF, from 129 Xe+ nat Sn, 32-50 AMeV), event by 
event properties at freeze-out were restored and used to build constrained caloric curves. The 
definition of pressure in the microcanonical ensemble is presented in section 2. Then, in section 
3, simulations to recover freeze-out properties of multifragmentation events |1U| [TT] are briefly 
described and constrained caloric curves deduced are discussed. Section 4 is dedicated to the 
description of the thermometer finally chosen and to comparisons of results relative to a caloric 
curve without any constraints. Constrained caloric curves with temperatures including quantum 
fluctuations are presented in section 5. Section 6 is devoted to a discussion of the various results. 
Conclusions are given in section 7. 



2. Pressure in the microcanonical ensemble 

Let us consider a gas of weakly interacting fragments (i.e. they interact only by Coulomb and 
excluded volume), which corresponds to the freeze-out configuration. Within a microcanonical 
ensemble, the statistical width of a configuration C, defined by the mass, charge and internal 
excitation energy of each of the constituting Mc fragments, writes 
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where I is the moment of inertia, K is the thermal kinetic energy, V is the freeze-out volume 
and xV M ° stands for the free volume or, equivalently, accounts for inter- fragment interaction 
in the hard-core idealization. 

The microcanonical equations of state are 
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Taking now into account that S = In Z = In J2c w c and that dW c /dV = (M c /V) W c , it 
comes out that 
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The microcanonical temperature is also easily deduced from its statistical definition |12j : 

= < (3/2M c -5/2)/K y- 1 . (4) 
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where K is the total thermal kinetic energy of the system at freeze-out. 

As Mo, the total multiplicity at freeze-out, is large, the pressure P can be well approximated 

by 

p _ 2 < g > 
3 V ' 



Z=6 Ebeam=32AMeV 




Z=6 Ebeam=39AMeV 

£200 
= 150 



Z=6 Ebeam=45AMeV 



Z=6 Ebeam=50AMeV 



10 20 » 

'frag 

Z=11 Ebeam=32AMeV 




10 20 



V frag? mlVl7ni ' 

Z=11 Ebeam=39AMeV 



3Q 40, 50 

1 



10 20 30 40 50 

V fragT mm/nS T 

Z=11 Ebeam=45AMeV 



10 20 30 40, 50 

V frag<. mm/nS ) 

Ebeam=50AMeV 




10 2 v° f -Tmm/nsT 

frag v ' 

Z=18 Ebeam=32AMeV 



=18 Ebeam=39AMeV 



frag < 

18 Ebeam=45AMeV 




Z=18 Ebeam=50AMeV 
in 15f 




10 

Z=27 Ebeam=39AMeV 
in 15 



10 20 

V frag 

Z=27 Ebeam=45AMeV 



v, ?mni/niy 



10 20 30 40, 50 

V frag( mm/nS )' 

Z=27 Ebeam=50AMeV 



c 
O 

o 



10 20 30 40, 50 

V frag( mm/nS T 







f 


: : : 


1 




m. 


10 2( 


> 30 40, 50 



c 

o 
O 



frag 



10 20 30 40, 5l 
V fragT mm/nS 



c 
O 

o 



J 



mi 



10 20 

v. 



30 40, 50 
', rag (mm/nsT 



Figure 1. Comparison between the experimental velocity spectra (full points) of fragments 
of a given charge and the simulated ones (histograms). Each row refers to a different fragment 
charge: starting from the top Z=6, Z=ll, Z=18 and Z=27. Each column refers to a different 
beam energy: starting from the left 32, 39, 45 and 50 AMeV. From 



3. Event by event freeze-out properties 

Starting from all the available asymptotic experimental information (charged particle energy 
spectra, average and standard deviation of fragment velocity spectra and calorimetry) of selected 
QF sources produced in central 129 Xe+ na *Sn collisions which undergo multifragmentation, a 
simulation was performed to reconstruct freeze-out properties event by event |10| 111) . The 
method requires data with a very high degree of completeness (total detected charge >93% of 
the total charge of the system), which is crucial for a good estimate of Coulomb energy. QF 
sources are reconstructed, event by event, in the reaction centre of mass, from all the fragments 
and twice the charged particles emitted in the range 60 — 120° in order to exclude a major part of 
pre-equilibrium emission. Dressed excited fragments and particles at freeze-out are described by 
spheres at normal density. Then the excited fragments subsequently deexcite while flying away. 
Four free parameters are used to recover the data at each incident energy: the percentage of 




Figure 2. Comparison betwwen the experimental (full points) and simulated (histograms) 
reduced relative velocity correlation functions for fragments. Each panel refers to a different 
beam energy. From 



measured particles which were evaporated from primary fragments, the collective radial energy, 
a minimum distance between the surface of products at freeze-out and a limiting temperature for 
fragments. The limiting temperature, related to the vanishing of level density for fragments |13j . 
was mandatory to reproduce the observed widths of fragment velocities. Indeed, Coulomb 
repulsion plus collective energy plus thermal kinetic energy (directed at random) plus spreading 
due to fragment decays are responsible for about 60-70% of the observed widths. By introducing 
a limiting temperature for fragments, thermal kinetic increases, due to energy conservation, 
which produces the missing percentage for the widths of final velocity distributions. The 
agreement between experimental and simulated velocity /energy spectra for fragments and for the 
different beam energies is quite remarkable (see figure [1]). Relative velocities between fragment 
pairs were also compared through reduced relative velocity correlation functions |144 [15] (see 
figure[2|). Again a good agreement is obtained between experimental data and simulations, which 
indicates that the retained method (freeze-out topology built up at random) and parameters 
are sufficiently relevant to correctly describe the freeze-out configurations, including volumes. 
Finally one can also note that the agreement between experimental and simulated energy spectra 
for protons and alpha-particles (see an example in figure [3]) is not very good at high energy. This 
comes from the fact that we have chosen (to limit the number of parameters of the simulation) 
a single value for the percentage of all measured particles which were evaporated from primary 
fragments. We shall see in the following section how to correct for temperature measurements 
derived from protons. 

From the simulations it is then possible to recover, event by event, the different quantities 
needed to build constrained caloric curves, namely the thermal excitation energy of QF hot 
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Figure 3. Centre of mass energy spectra for protons, deuterons, tritons, 3 He, 4 He and 7 Li 
for the reaction at 39 AMeV. Open points are the simulated data, while full points are the 
experimental ones. From 



nuclei, E* (total excitation minus collective energy) with an estimated systematic error of around 
1 AMeV, the kinetic temperature T^ in at freeze-out, the freeze-out volume V (see envelopes of 
figure 8 from and the total thermal kinetic energy at freeze-out K. In simulations, Maxwell- 
Boltzmann distributions are used to reproduce the thermal kinetic properties at freeze-out and 
consequently the deduced temperatures, Xfcj n , are classical. 

Constrained caloric curves, which correspond to correlated values of E* and T^ n have been 
derived for QF hot nuclei with Z restricted to the range 80-100, which corresponds to the A 
domain 194-238, in order to reduce any possible effect of mass variation on caloric curves [7] 
(see figure 0|). Curves for internal fragment temperatures are also shown in the figure. For 
different average freeze-out volumes expressed as a function of Vo, the volume of the QF nuclei at 
normal density, a monotonous behaviour of caloric curves is observed as theoretically expected. 
The caloric curves when pressure has been constrained exhibit a backbending and moreover 
the qualitative evolution of curves with increasing pressure exactly corresponds to what is 
theoretically predicted with a microcanonical lattice gas model [8]. However by extrapolating at 
higher pressures, one could infer a critical temperature around 20 MeV. Such a value is within 
the range calculated for infinite nuclear matter whereas a lower value is expected for nuclei 
in relation with surface and Coulomb effects pQ. At that point, one may wonder if Tki n is a 
relevant thermometer. We shall see, in what follows, that the final choice was to use, for protons 
thermally emitted at freeze-out, a new thermometer recently proposed for which quantum effects 
can be included. 
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Figure 4. Caloric curves (kinetic temperature versus thermal excitation energy) constrained at 
average volumes (left) and for selected ranges of pressure (right) and the corresponding internal 
temperatures for fragments. Error bars correspond to statistical errors. From [16j . 



4. Temperatures from proton transverse momentum fluctuations 

Very recently a new method for measuring the temperature of hot nuclei was proposed [17U18j. It 
is based on momentum fluctuations of emitted particles like protons in the centre of mass frame of 
the fragmenting nuclei. On the classical side, assuming a Maxwell-Boltzmann distribution of the 
momentum yields, the temperature T is deduced from the quadrupole momentum fluctuations 
defined in a direction transverse to the beam axis: 
v 2 = < Ql y > - < Qxy > 2 = 4m 2 T 2 

with Q xy = P 2 - P 2 ; m and P are the mass and linear momentum of emitted particles. 
Taking into account the quantum nature of particles, a correction Fqc related to a Fermi- 
Dirac distribution was also proposed |18[ [TU] . 
In that case a 2 = 4m 2 T 2 F QC where F QC = 0.2(T/e / )~ L71 + 1; 

ef = 36 (p/po) 2 / 3 is the Fermi energy of nuclear matter at density p and po corresponds to 
normal density. 

Before using the thermometer with protons to build constrained caloric curves, it was 
important to verify several things. With the classical simulation (freeze-out and asymptotic 
proton momenta) , it is possible to test the agreement with the proposed classical thermometer. 
Moreover the effects of secondary decays on temperature measurements can be estimated. 
Figure [5] shows different caloric curves without constraints. Note that the selection in Z and 
A of hot nuclei is the same as in the previous section. It was verified that, within statistical 
error bars, at a given thermal excitation energy transverse momentum fluctation values are the 
same for our selection or by selecting only a single (A and Z) hot nucleus. Open diamonds refer 
to classical temperatures calculated from momentum fluctuations for protons thermally emitted 
at freeze-out. Within the statistical error bars they perfectly superimposed on T^i n values 
(see figure 2 of [IB])) which just verifies that Maxwell-Boltzmann distribution was correctly 



implemented in the simulation. Full squares correspond to classical temperatures calculated 
from momentum fluctuations for protons after the secondary decay stage. We note that the 
caloric curve is distorted, which means that it is not possible to use experimental data from 
protons to measure temperatures. Moreover, in this case, quantum corrections for temperatures 
can not be made because protons are emitted at different stages of deexcitation with different 
Fermi energy values. In figure [5] classical temperatures calculated from experimental proton 
data are also shown (full points). As for temperatures calculated from asymptotic proton data 
of simulations, a monotonous behaviour of the caloric curve is observed. One also note the 
differences between the two sets of temperature values, which are related to the fact that, as 
indicated in the previous section, simulations do not describe experimental proton energy spectra 
very well. Those temperature differences will be used to correct classical temperatures derived 
from simulated protons at freeze-out. 
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Figure 5. Caloric curves (classical tem- 
perature from proton transverse momen- 
tum fluctuations versus thermal excitation 
energy) for protons (simulation) thermally 
emitted at freeze-out (open diamonds), for 
protons (simulation) after the secondary 
decay stage (full squares), and from pro- 
tons experimentally measured (full points). 
Error bars correspond to statistical errors 



Figure 6. Caloric curves: classical 
temperature (open diamonds)/ quantum 
corrected temperature (full squares) from 
proton transverse momentum fluctuations 
versus thermal excitation energy. Protons 
(simulation) are thermally emitted at 
freeze-out. Error bars include statistical 
and systematic errors. 



It finally appears that the only way to extract temperatures from proton transverse 
momentum fluctuations taking into account quantum effects is to use protons thermally emitted 
at freeze-out. In that case classical temperature values from simulations must be extracted 
and corrected and then, quantum corrections applied, which needs Fermi energy values. Those 
values can be estimated from semi-classical calculations (Xn+Sn at 32 AMeV and Sn+Sn at 50 
AMeV) |20[ I21j: protons thermally emitted at freeze-out at time around 100-120 fm/c after the 
beginning of collisions come from a low density uniform source. For the two incident energies 



low densities around p ~ 0.4po are calculated which corresponds to e/ ~ 20 MeV. We have 
introduced a systematic error of ± O.lpo f° r the calculation of e/ and consequently a systematic 
error for "quantum" temperatures of ± 0.6-0.5 MeV on the considered thermal excitation energy 
range. Figure [6] shows the final caloric curve with temperatures from quantum fluctuations (full 
squares). It exhibits a plateau around a temperature of 10-11 MeV on the E* range 5-10 AMeV. 
For comparison the caloric curve with classical temperatures derived from the simulation is also 
shown (open diamonds). 



5. Constrained caloric curves with quantum temperatures 

Constrained caloric curves, which correspond to correlated values of E* and quantum corrected 
temperatures have been determined. Accordingly E* values, which were initially derived from 
experimental calorimetry including estimated correction for neutrons (see [10]). have been 
corrected a posteriori using quantum temperatures at freeze-out. Pressure values were also 
corrected using quantum temperatures. In figure [7] (left) we have constructed caloric curves for 
two different average freeze-out volumes corresponding to the ranges 3.0-4.0Vb and 5.0-6.0Vo 
where Vo correspond to the volume of the QF nuclei at normal density. Again as theoretically 
expected a monotonous behaviour of caloric curves is observed. Figure [7] (right) shows the 
caloric curves when pressure has been constrained within two domains: 1.3-4.5 and 4.5-7.9 10~ 2 
MeVfm~ 3 . Backbending are clearly seen especially for the lower pressure range. For higher 
pressures the backbending of the caloric curve is largely reduced and one can estimate that the 
critical temperature is around 12-13 MeV for the selected hot nuclei. 
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Figure 7. Caloric curves (quantum corrected temperature versus thermal excitation energy) 
constrained at average volumes (left) and for selected ranges of pressure (right). Error bars 
include statistical and systematic errors. 
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Figure 8. Caloric curves (classical temperature derived from experimental proton data versus 
thermal excitation energy) constrained at average volumes (left) and for selected ranges of 
pressure (right). Error bars correspond to statistical errors. 



6. Discussion 

First of all, as far as internal temperature for fragments are concerned (see figure H]), one 
can observe that they perfectly agree, for our A range, with those calculated with the well 
known "He/Li thermometer" used in [5j[7] keeping the proposed prefactor 16. Note that similar 
temperature values are also obtained from our data (see figure 2 from [E]). This is a strong 
indication that the He/Li thermometer mainly reflects the internal temperature of fragments 
in relation with the strong secondary emission of alpha-particles and Li isotopes. As indicated 
in [7j those plateau temperatures can be interpreted as representing the limiting temperatures 
resulting from Coulomb instabilities of heated nuclei which has long been predicted [22] , Indeed, 
for thermally equilibrated QF hot nuclei one expects internal temperature for fragments equal to 
the limiting temperature of the fragmenting system. As a direct consequence internal fragment 
temperatures must reflect the evolution of limiting temperature with A of hot nuclei, which is 
indeed experimentally observed [7|. 

For a finite piece of nuclear matter like a hot nucleus the micro canonical ensemble is the 
most relevant ensemble and the kinetic/microcanonical temperature is the relevant parameter to 
build caloric curves and deduce information on a possible phase transition. On the other hand, 
one may also recall that dispersions in fragment transverse momentum spectra generated in 
projectile fragmentation were succesfully explained by adding the Fermi momenta of individual 
nucleons in fragments [23], which is a consequence of the quantum nature of nucleons. For 
multifragmentation reactions a similar approach was first proposed in [23] to qualitatively 
explain the different temperatures obtained from various observables (kinetic energies, level 
populations, isotope ratios). Thus, the use of the thermometer recently proposed based on 
quantum fluctuations [18] was a good opportunity to better investigate the behaviour of caloric 
curves especially under constraints in volume and pressure. The behaviour predicted by 



theoretical works for a first order phase transition is observed and confirms results obtained 
for other signals: negative microcanonical heat capacity and bimodality of an order parameter. 

The last point that we want to discuss concerns information which can be directly deduced 
from experimental data like proton transverse momentum fluctuations. In that case as previously 
indicated Fermi energies can not be estimated and only classical temperatures can be calculated. 
From simulations we have seen previously that secondary decays distorted caloric curves. Results 
displayed in figure [H] fully confirm this observation: constrained caloric curves built with classical 
temperatures derived from experimental proton momenta also exhibit a monotonous behaviour 
and do not show any dependence upon average volumes or pressures, which prevent any direct 
information. 



7. Conclusions 

Several caloric curves have been derived for hot nuclei from quasi-fused systems using a new 
thermometer based on proton transverse momentum fluctuations including quantum effects. 
The unconstrained caloric curve exhibits a plateau at a temperature of around 10-11 MeV on 
the thermal excitation energy range 5-10 AMeV. For constrained caloric curves (volume and 
pressure) we observe what is expected for a first order phase transition for finite systems in 
the microcanonical ensemble, namely a monotonous behaviour at constant average volumes and 
backbending for constrained pressures. After the observation of negative microcanonical heat 
capacity and bimodality the behaviour of caloric curves is the ultimate signature of a first order 
phase transition for hot nuclei. 

The exit, at high excitation, of the spinodal region and of the coexistence region around 
respectively 8 and 10 AMeV are in good agreement (within error bars) with spinodal [25] and 
bimodality [I] signals. Note also that 10 AMeV corresponds to first experimental indications 
for the onset of vaporization [26] 1271 1281 129] . 

Finally, one can infer from the present studies that the temperatures obtained with the He/Li 
thermometer, which was largely used in the past, seem mainly reflect the internal temperatures 
of fragments in the excitation energy range 5-10 AMeV. 
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